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Abstract: 

The dispersion relation for the inverse hyperbolic potential is calculated in the 
classical limit. This is shown for both the low amplitude phonon branch and the 
high amplitude soliton branch. It is shown these results qualitatively follow that 
previously found for the inverse squared potential where explicit analytic 
solutions are known. 
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Introduction 


In one spatial dimension, a class of integrable many-body systems is the 
Calogero-Sutherland-Moser systems. They consist of many identical non- 
relativistic particles interacting through two-body potentials of the inverse 
square type and various extensions such as the inverse sine squared, the inverse 
hyperbolic squared, and the inverse squared Jacobian elliptic function^' 2 ' 3 ). The 
significance of these models is they represent a class of integrable many-body 
systems. One test for integrability is to use the Lax method, which is the method 
first applied by Calogero, et. al to the inverse squared potential®. Using this 
method one attempts to show that for certain potentials one can find two 
Hermitean N x N matrices L and A that follow the Lax equation dL/dt=i ( AL-LA ). 
With this L evolves as a transformation generated by A and det[ L-/.1] is a 
constant of motion. 

In this paper we look at the inverse hyperbolic squared potential and 
derive the dispersion relation in the classical limit. This is useful because it is one 
method one can compare the quantum mechanical system with the purely 
classical many-body soliton system, which may also be derived using the 
classical equations of motion. This comparison was done in a previous paper for 
the inverse squared potential® . However, the dispersion relation for this system 
is a bit more problematic as no closed form analytic solution is known. 

Dispersion relation for the inverse hyperbolic potential 

We begin with the many body Hamiltonian 
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We make a change of variable, x ox, so that 1/a is the unit of length, and so 
the Hamiltonian becomes 
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where energies are measured in units of h 2 a 2 / m, and X(X-l)=gm/2/z 2 . Thus, we 
see that the classical limit h-> 0 corresponds to A — > oo ; this is the limit we 
investigate in this paper. This quantum system was originally solved by 
Sutherland^' 2 ' 3 ), using the asymptotic Bethe ansatz^TT), an( j a proof of 
integrability due to Calogero, et. al.W 

As is well knowndTT), the solution by the asymptotic Bethe ansatz for the 
low energy properties makes use of two quantities p(k) and e(k), which are 
solutions of the following integral equations: 

P (k) + -L r B e ' (k - k’ ) P (k' )dk* = 0) 

Z7lJ-o 2n 




and 



The limit to the integral B is related to the density d=N/L through 


(4) 


j^ B p(k)dk=d, 


(5) 


while the chemical potential ii is determined by the requirement that e(±B)=0. 
The kernal 0’(k) of the integral equation is the derivative of the two-body phase 
shift 0(k), which for our (reduced) Hamiltonian is 


0(k) = ilog 


T(1 + ik / 2)T(X - ik/ 2) ~ 
T(1 — ik / 2)T(X, + ik/ 2) ' 


(6) 
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The low-lying or zero temperature physical properties are determined from p(k) 
and e(k) as follows. First, the ground state energy is calculated as 

E/L =i/-B p(k)k2dk - (7) 

From this, the zero temperature equation of state may be determined. Second, 
the dispersion relation for low-lying excitations may be determined 
parametrically, with the energy given by I s(k) I , and the group velocity v(k) by 
v(k)=E'(k)/2jtp(k). (We choose the group velocity rather than the momentum to 
make the classical limit easier.) The dispersion relation has two branches: 

I k I <B, called the hole or phonon branch; and I k I >B, called the particle or 
soliton branch. This is discussed by Sutherland*^ 1 ), where the quantum Toda 
lattice is first solved, and the classical limit taken to recover the original results of 
Toda( 7 ). 

In the classical limit A -» we find B -» oo, and so we rescale k by k=2Ax. 
Then, the kernal of the integral equation becomes 

0’(k)-ilog[l + l/x 2 ]. (8) 

Letting B=2Ab, and keeping terms of leading order in A, we find that for I x I <b, p 
and e obey the integral equations 

J^ b dy log[l + 1 / (x - y) 2 Ap(2 Ay) = 1 , (9) 

and 

J" b dylog[l+ l/(x -y) 2 A _1 e(2Ay) = 4itx 2 -2tc p/ A 2 . (10) 
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We see that the quantities that approach finite limits in the classical limit are 
Xp(2Xx)=p(x) and X' 1 e(2Xx)=s(x) / for I x I <b. Outside the integration region, when 
I x I >b, we use the full integral equations, so that 

p(2Xx) = ^--^-J^dylog[l+l/(x-y) 2 p(y) (11) 

and 

X~ 2 e(2Xx) = 2k 2 - ^ - ^-J^dylog^l + l/(x - y) 2 e(y) . (12) 

Thus, on the soliton branch, we define the quantities which approach finite 
classical limits as p(2Xx)=p(x) and X' 2 e(2Xx)=e(x), for I x I >b. 

Looking at the group velocity, v(k)=s'(k)/2jrp(k), and using 
e'(2Xx)=E'( x )/2, | x I <b, or X"1 e'(2Xx)=e'( x )/ 2, lxl>b, we see that X" 

1 v(2Xx)=£ , ( x )/4jrp(x)=v(x), for all x. Thus, the dispersion relation is given in the 
classical limit by velocity Xv(x) and energy X I e(x) I , I x I <b, or X 2 I e(x) I , I x I <b. 
This gives an energy proportional to h for the phonon branch, as is to be 
expected. 

We thus see that the dispersion relation in the classical limit depends on 
the solution of the two integral equations 

dy log[l + 1 / (x - y) 2 p(y) = 1 (13) 

and 

J ^ dylog 1+ l/(x -y) 2 e(y) = 4jtx 2 -2jtp. (14) 

In the second equation, u represents the old u/X 2 . The phonon branch of 
the dispersion curve can be found by a straightforward harmonic approximation 
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to the equations of motion. By familiar methods, this gives for the frequency 
oo(k) of an oscillation of wavenumber k. 


oo 2 (k) = 4A 2 2 [1 -cos(kj/ d)] 
j=l 


3 + 2 sinh 2 (j / d) 
sinh 4 (j/ d) 


(15) 


Here, d is again the density N/L. The group velocity is given by the usual 
expression v(k) = dio(k) / 3k, so the dispersion relation for the phonon branch w(v) 
at density d is thus determined parametrically. This must coincide with the 
previous expression when I v I <v s = v(0). 

To check this correspondence, we make what are apparently convergent 
expansions for both p(x) and e(x). 


p(*) = 


cl 


In-^b 2 - . 


1+2 PjT 2J (x/b) 


7-1 


(16) 


and 
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-x 2 S e jU 2 j (x/b). 

j-0 


(17) 


Here Tj(x) and Uj(x) are Chebyshev polynomials of the first and second kind 
respectively. Using these expansions, which appear rapidly convergent except 
for the limit b -> oo, corresponding to the inverse square potential, which can 
independently evaluated^ 8 ' 9 ), one finds the dispersion relations shown in figure 
(1) for selected values of the density. This agrees numerically for the phonon 
branch with the previous harmonic approximation. 
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Figure 1 . We show the dispersion relations for selected values of the density d = 
3.058..., 1.729..., 1.039..., and 0.673.... These values interpolate between the 
inverse square interaction and the Toda lattice. To the left we show the phonon 
dispersion relation by plotting the frequency as a function of the group velocity, 
and to the right we show the soliton dispersion relation by plotting the energy of 
the soliton as a function of the soliton velocity. 

Discussion 

The dispersion relation qualitatively follows that found previously for the 
inverse squared potential in the classical limit! 5 ). In that case simple analytic 
results were found. For each density shown, the point where the curve for both 
the soliton branch and the phonon branch is zero represent where the velocity is 
the speed of sound for that density. The phonon branch represents simple low 
amplitude sound waves whereas the soliton branch represents high amplitude 
non-linear waves. We have found these results by taking the classical limit of a 
quantum mechanical system. It would be interesting to see how this compares to 
the direct solution of the classical equations of motion as we have done 
previously for the inverse squared potential! 5 ). We plan to investigate this in a 
future work. 
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